ASYMPTOTICS OF MAXWELL TIME IN THE PLATE-BALL PROBLEM 
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Abstract. The problem on rolling of a sphere on a plane without slipping or twisting is considered. 
One should roll the sphere from one contact configuration to another so that the length of the curve 
traced by the contact point in the plane was the shortest possible. Asymptotics of Maxwell time for 
" rolling of the sphere along small amplitude sinusoids is studied. Two-sided estimate for this asymptotics 

I . is obtained. 

o; 
<n: 

1. Introduction 

For the problem on rolling of a sphere on a plane without slipping or twisting, an optimal control 
\Q \ problem is studied. State of the system is described by the contact point of the sphere and the plane 
and orientation of the sphere in three-dimensional space. One should roll the sphere from one contact 
i—i. configuration to another so that the length of the curve traced by the contact point in the plane was the 
U: shortest possible. The problem has application in robotics: rotation of a solid body in robot's hand. 

In this work we obtain two-sided estimate of Maxwell time in the plate-ball problem in asymptotic 
c| ■ case. 

The problem was stated in [1] by Hammersley. Then Arthur and Walsh [2] proved integrability 
of Hamiltonian system of PMP in elliptic functions. Jurdjevic in [3, 4] showed that projections of 
'— extremal curves (x(t),y(t)) are Euler elasticae (see the papers [5, 6])). He gave a description of 
I . different qualitative types of extremal trajectories, and obtained differential equations for evolution 
of Euler angles along extremal trajectories. Explicit formulas for the extremals were obtained in the 
paper [7]. 

. Optimality of extremals is still an open problem nowadays. Short arcs of extremal trajectories are 
optimal but long arcs, in general, are not optimal. A point at which an extremal trajectory loses global 
optimality is called a cut point. A cut point is a conjugate point or a Maxwell point. A Maxwell point 
is a point in the state space, where an extremal trajectory crosses another one with the same value of 
cost functional. Yu. Sachkov began to study cut points in the plate-ball problem (see the paper [8]). He 
found continuous and discrete symmetries of the exponential mapping. Then he obtained equations, 
which define Maxwell points as fixed points of the discrete symmetries, and formulated necessary 
optimality conditions in terms of Maxwell time (see Theorems A, B in Section 2). But the problem 
of optimality of extremal trajectories is steel open because the equations, which define the Maxwell 
points, are not solved. 

The papers [7, 9] present the asymptotics of extremal trajectories in a neighborhood of the stable 
equilibrium of a mathematical pendulum (see (6)), which appears in the adjoint subsystem of the 
Hamiltonian system of the Maximum Principle. In this case, the extremal curves on the plane are 
close to sinusoids of small amplitude. 

This work continues to study the problem of optimality of extremal trajectories. It studies the 
problem in an asymptotic case, where the formulas defining the extremal trajectories and the Maxwell 
points are expressed via trigonometric functions. They are simpler than in the general case where 
the formulas are expressed in elliptic functions. We study the behavior of Maxwell points MAX 1 
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and MAX 2 in this case and obtain two-sided estimates for the first Maxwell times t\ and t 2 , which 
correspond to the fixed points of the discrete symmetries e 1 and e 2 of the exponential mapping. 

2. Statement of the problem and known results 

In this section we formulate the plate-ball problem and recall some known results. Let (x,y) G E 2 
be the contact point of the sphere and the plane. By q = (90, 9i, 92, 93) G S 3 denote the unit quaternion 
(see the paper [10]) representing the rotation of three-dimensional space, which translates the current 
orientation of the sphere to the initial orientation. The problem on optimal rolling of a unit sphere 



on a plane is stated as follows: 

Q = u 1 X 1 (Q)+u 2 X 2 (Q), (1) 

Xi(Q) = (1,0, q 2 ,q 3 , -q ,-qif, X 2 (Q) = (0, 1, - Ql , q Q ,q 3 , -q 2 ) T , (2) 

Q = (x,y, go,9i,92, 9s) G M = E 2 x S 3 , u = (m, u 2 ) G E 2 , (3) 

Q(0) = Q = (0,0,1,0,0,0), Q(h) = Qi, (4) 

Z = J \J u\ + u| — > min . (5) 
o 



Admissible controls are measurable and essentially bounded. Admissible trajectories are Lipschitz. 
Problem (l)-(5) is a left-invariant sub-Riemannian problem on the Lie group M = E 2 x S 3 . The control 
system is completely controllable by the Rashevsky-Chow theorem (see the paper [11]). Existence 
of optimal controls follows from Filippov's theorem (see [11]). The Maximum Principle is applied to 
study the optimal controls. In the abnormal case, the sphere rolls along a straight line in the plane 
(x, y). In the normal case, the subsystem of the Hamiltonian system for the adjoint variables (9, c, r, a) 
satisfies the equations of a mathematical pendulum as follows: 

9 = c, c=— rsin#, a = r = 0. (6) 

Projections of extremal trajectories to the plane (x, y) are Euler elasticae, i.e. stationary configurations 
of an elastic rod on a plane with fixed end points and fixed tangents at these points (see the paper [3]). 

In the paper [8] the author describes continuous and discrete symmetries of the exponential mapping 
in the plate-ball problem 

Exp : (X,t) (->■ Q t , (X,t) E N = C x E + , Q t G M = E 2 x SO(3), 

C = {A G T«5 M | H{\) = 1/2} = {(0,c, a, r) \ 9 G S 1 , c G E, r ^ 0, a G 

The continuous symmetries \ (3 G S 1 } are rotations by the angle /3 in the plane (x,y). The discrete 
symmetries e 1 , e 2 , e 3 are reflections of the trajectories of pendulum (6) in the coordinate axes {c = 0}, 
{9 = 0}, and in the origin (6>, c) = (0,0) respectively. The action of the symmetries in the preimage 
and in the image M of the exponential mapping is defined in [8]. Also, there is a description of the 
Maxwell sets corresponding to the symmetries e l , i = 1, 2, 3: 

MAX* = {(A,t) G N | 3/3 G S 1 : (\,t) =e i o^(\,t), Exp(A,s) ^ Exp(A,s), Exp(A,i) = Exp(A,t)}. 
In particular, there is proved the following theorem for the symmetry e 1 . 

Theorem A ([8]). Suppose t > and Q s = (x s ,y s ,R s ) = Exp(A,s) is an extremal trajectory such 
that the elastica {(x s ,y s ) \ s G [0,t]} is nondegenerate, is not centered at an inflection point, and 
satisfies the following equation: 

? 3 (*)=0. (7) 
Then (A,t) G MAX 1 . Therefore the trajectory Q s , s G [0,t] is not optimal for any t > t. 
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There is a similar theorem for the symmetry e 2 . 

Theorem B ([8]). Suppose t > and Q s = (x s ,y s ,R s ) = Exp(A,s) is an extremal trajectory such 
that the elastica {(x s ,y s ) \ s G [0, t]} is nondegenerate, is not centered at a vertex, and satisfies the 
following equation: 

(x qi +yq 2 )(t)=0. (8) 
Then (A,i) G MAX 2 . Therefore the trajectory Q s , s G [0, t] is not optimal for any t > t. 

In Sections 4-7 we study asymptotics of the first Maxwell times ti and t 2 corresponding to the 
Maxwell points (A,t) G MAX 1 and (A,t) G MAX 2 , near stable equilibrium of mathematical pendu- 
lum (6). Note that there is a similar theorem for the the symmetry e 3 in [8], but we do not consider 
the Maxwell set MAX 3 since it has lesser dimension than MAX 1 and MAX 2 . 



3. Asymptotics of extremal trajectories and 
limit behavior of Maxwell sets MAX 1 and MAX 2 

We study optimality of extremal trajectories in the plate-ball problem. Due to complexity of 
parametric equations of these trajectories it is very difficult to study this problem in the general case. 
Therefore we study the asymptotic case, corresponding to small oscillations of pendulum (6). In this 
case, the sphere rolls along the curves close to small-amplitude sinusoids. 

We consider the Hamiltonian system of the Maximum Principle near the stable equilibrium 9 = c = 
of mathematical pendulum (6). The asymptotics as p\ = 6q + c\ — > of the solutions x(t), y(t), qo(t), 
Qi(t)i Q2(t), qs(t) of this system is presented in [9]. It is well known that mathematical pendulum (6) 
is a harmonic oscillator in the asymptotic case. The Hamiltonian system is simplified by change of 
variables corresponding to the symmetry of the system "rotation by angle a", which is defined as 
follows: 

(t, 9, c, a, r, x, y, ui , u 2 , qo, qi, q 2 , qs) -> (s, 0, d, a, m, x, y, ui, u 2 , %, qi,q 2 , q 3 ), 
s = mt, d = c/m, m = ^fr, 

x ^ = A(a) ( l), I": )=A(a)( ? ), l* = «" (9) 



cos a - 


sin a 


sin a cos a 


I)' \ 






[Q3 : 



y ) \y ) v 92 / v 92 / [g 3 = q 3 

Asymptotics for the elastica x, y and the components qi of the quaternion q are the following: 

x(s) =- + 0( p 2 ), y(s) = l(0 osins + d (l- cos s)) + 0(p 2 ) (10) 
mm 

9o(s) = cos ^-+O(p 2 ), q 2 (s) = -sm^- + O(p 2 ), (11) 
= 2{m l_ l) (m cos J- sin s - (1 + cos s) sin ^-)9 0+ 

+ 2(mi-l) ( m( - 1 ~ C ° S ^ C ° S 2m" ~ Sin S Sin 2m" ^° + ^ 
^ {S) = 2(mi-l) i{ - 1 + C ° S S) C ° S 2m" + m Sin S Sin 2m"^ 0+ 

+ 2(m^T) (sin S COS 2m" " m(1 + C ° S S) + (13) 
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These expansions have a removable singularity at m = 1. Up to O(pq) the curve (x,y) is a sinusoid 
of small amplitude Complexity of the formulas for the asymptotics of the quaternion is related 
to different frequencies of trigonometric functions. This implies existence of "resonance" instants. 
At such instants plots of the asymptotics of the first Maxwell times t\ and £2, corresponding to the 
symmetries e\ and £2 respectively, have vertical tangent lines, see below. In the following sections we 
use (10)-(13) to study asymptotics of Maxwell points as po — > 0. 

To examine optimality of extremal trajectories in asymptotic case we study the first Maxwell times 
ti and ti- They are first values of time, when extremal trajectory arrives at the Maxwell sets MAX 1 
and MAX 2 respectively. 

By Theorem A the equation c/3 = implies existence of a Maxwell point for the elasticae that are 
nondegenerate and not centered at an inflection point. In view of (9)-(13), the leading term of the 
function qs(t) has a root if there holds the following equality: 

4 cosp -0 O ship p . . P> „ s 

t. (cos — snip — mcospsin — ) = 0, p=-. 

m 2 — 1 m m 2 

Roots of the factor do cosp — 0o ship has a simple geometric meaning for the sinusoid (x°(s), y°(s)) = 
(s/m, (0osins + do(l — coss))/m) (the leading term of the asymptotics for the elastica (x(s),y(s))), 
and therefore for the sinusoid (x° (s) , y° (s)) = (cos a x°(s) + sin a y (s), — s'max°(s) + cos a y°(s)) (the 
leading term of the asymptotics for the elastica (x(s),y(s)) as po — > 0). It can easily be checked that the 
sinusoid {(x° (a) , y° (a)) \ a € [0, s]} is centered at the inflection point if and only if do cosp — #0 sin p = 
0. 

Thus we study roots of the factor (cos — sinp — mcospsin ^)/(m 2 — 1). This factor has an unre- 
movable singularity if m = 0, and it does not vanish for any p 7^ if m = 1. Therefore we consider 
the function 

p p 
Oi(p,m) = cos — sinp — mcospsin — , m € (0, 1) U (1, +00), p > 0, (14) 
m m 

and study its minimal positive root 

Pi(m) = min{p > 0\gi(p,m) = 0}. (15) 

In Section 4 we study the function pi(m). We obtain two-sided estimate of pi(m) and prove that 
Pi(m) is monotone and differentiable for m > 0, m 7^ 1. 

Similarly, we examine the first Maxwell time ti (the first instant of time when an extremal trajectory 
rich the Maxwell set MAX 2 ). By Theorem B the equation xq\ + yq2 = implies existence of a Maxwell 
point for the elasticae that are nondegenerate and not centered at the top. In view of (9)-(12), the 
function xq\ + yq2 is equal to zero up to O(pq) if there holds the following equality 

do sinp + O cosp . . s 

7 — 2 r\ 92{P,m) = 0, P=-, 

m(m z — 1) 2 

where 

92(p, m) = mpcos (—) sinp — (pcosp + (m 2 — 1) sinp) sin(— ). (16) 

If an elastica is not centered in the vertex, then the factor do sinp + 0ocosp is bounded away from 
zero. Therefore, we are interested in a minimal positive root of the function g2(p,m) = 0: 

P2(m) = min{p > 0\g2(p, m) = 0}. (17) 

In Section 5 we study the function P2(m). We obtain two-sided estimate of P2{m) and prove that 
P2(m) is a continuous function for m > 0, m 7^ 1. 
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4. Study of the function p\(m) 

In this section we prove a two-sided estimate of a minimal positive root of the equation g\ (p, m) = 
(see (14)). Our aim is to find or estimate as accurately as possible for any m G (0, 1) U (l,+oo) the 
value of pi(m) defined in (15). The main results on this problem are summarized in Theorem 4.1. 

Theorem 4.1. For any m G (0, 1) U (1, +oo) the function g\(p, m) has a minimal positive root pi(m) 
satisfying the following properties: 

a) The function pi(m) is continuous and increasing for m G (0,1), coincides with at the 
points {m = pp^A; € N}; is continuous and decreasing for m G (1, +oo), coincides with at 
the points {m=^|l;eN}. 

b) The function pi{m) is continuously differentiable at all points of the interval (0,1), except for 
the set {m = jjh[\k G N}, where its derivative is equal to +oo; is continuously differentiable at 
all points of the interval (1, +oo), except for the set {m = ^-\k G N}, where its derivative is 
equal to -co. 

c) For m G (|, 1) the plot of p = p\(m) "wraps" the hyperbole p = as follows: 



PM) < T^Tn forme (f^, 



k 
k+1 



(18) 

i=m <Pi(m)forme (jfe,^) , k G N. 
For m G (1,3) the plot of p = p\{rn) "wraps" the hyperbole p = -^zi as follows: 

d) Estimates (18), (19) are supplemented by estimates from the other side. For\/k G N define 



ftiH = min (#^) ~k,pk + l- arcsinl^j ,m G (§g, ^ , 

h 2 (m) = min (^i±s - (2* + 1), fc + 1 - ^ ^ arcsin , ro € (^, f±±) . 

The following inequalities are valid for m G (|, l) : 

' ~ gAiW < PiM /or m G fe, ^ 

TTte following inequalities are valid for 1 < m < 3: 

< + M m G (f±?, ^±1 

rr - TwMh) < piM V m g (^±i. . fc g n. 

e) For < m < | i/ie function pi(m) has the following estimates: 

, 7rm 37rm 

maxjpm, } < pi(m) < — — , (22) 

1 — m 2 

where p = 4.493409 ... is the root of the equation tanx = x that lies in the range it < x < 4^. 
For m > 3 the function p\(m) has the following estimates: 

, irm , . . 37T 
max{p, — — } < Pl (m) < — . (23) 

777* X ^ 

f) At the point p = pi(m) the function p i-> gi(p, m) changes its sign. 



(20) 



(21) 
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Fig. 1. Plot of the function pi(m) with two-sided estimates 



Figure 1 shows the plot of p\{m) with two-sided estimates (dashed lines). 

Remark 4.1. If m G [ffejj, feXj]) then the values of s(m) = range over the interval [2k, 2k + 1], 
s (ff^xj = 2/c, s (fpt) =2k + l. Consequently, 

^ hi(m) < 0.3264, lim 7ti(m) = 0, lim /ii(m) = 0. 
If m G [fc^j; 2^+3 ]' ^ en ^ e values of s(m) range over the interval [2k + 1, 2k + 2], s ( ^+3 ) = 2/c + 2. 



h 2 (m) < 0.3244, lim /i 2 (m) = 0, lim h 2 {m) = 0. 



Remark 4.2. From (18)-(21) it follows that pi(m) tends to +00 asm-s-110. 

4.1. Estimation of pi(m) for < m < |. Let us prove that equation (14) has a root in the interval 
< p < -^rp < |. In this interval the function sinp does not vanish. The function sin vanishes 
only at p = nm but this point is not a root of equation (14), since gi(Trm,m) = —sm(irm) 7^ 0. 
Therefore, dividing the equation g\ (p, m) = by their product we obtain in the interval < p < 5 
equivalent equation 



cot I — I = m 



cotp 44> / ( — ) = f(p), where f(x) = xcot(x). (24) 
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Put y = ^. Then equation (24) has the form f(y) = f(my). The function / decreases over the 
intervals (0,7r), (tt, 2tt) and satisfies the following equalities: 

3tt 

lim f(y) = 1, lim f(y) = -co, lim /(y) = +00, /(— ) = 0. 

It follows from decrease of f(y) over the interval < y < tt that the equation f(y) = f(my) has no 
roots for < m < 1 in this interval. But a root exists in the interval tt < y < ^ since the difference 
f(y) — f(my) is continuous and tends to +oo as y — > tt + and is negative at ^ (because of condition 
m < | we have / (^jp) > 0). Denote this root by yo- Since f{my) < 1 for < y < ^, it is clear 
that yo(m) exceeds the number p defined by the equation f(p) = 1, -it < p < Consequently, 
p < yo(m) < t^. Thus for < m < | we obtain the inequality: 

. , 3-7rm 

pm< Pl (m) < , (25) 

where p = 4.493409 ... is the root of the equation tanx = x that lies in the interval tt < x < 

From (25) it follows immediately that for < m < | the root pi(m) lies in the interval < p < |, 
as it was claimed at the beginning of this section. The lower estimate in (25) needs to be improved 
as m — > | — 0. This improvement will be done below in Section 4.3. 

4.2. Reducing the problem to finding (estimation) of the minimal positive root of a 
simpler equation g(x,s) = 0, s > 1, g(x,s) = ssinx — sin(sx). It can easily be checked that 

gi(p, m) = (sinpcos £ — cosnsin £-) — (sinpcos £ + cospsin £■) . 

Zi V Tit 77% / A V 77% 771 / 

Since sinpcos £ ± cospsin £ = sin (p ± £), we see that the function g\ admits the following repre- 

lib lib \ lit/ 

sentation: 

9l (p, m) = 2^1 sin (p (1 - i)) - 2^1 sin (p (l + ±)) = ^ sin (p2S±l) _ i±™ sin (pl=p) . 
Let us introduce a parameter s = and a new variable x = P^—^- Now express gi(p,m): 

1 — m 1 + m 1 — m m — 1 

gi{p,m) = — - — sm(sx) — smx = — - — (sin(sx) - ssinx) = — - — g(x,s). 

Thus, we study the function 

x±(s) = min{x > 0\g(x,s) = 0} for s > 1. 

4.3. Absence of roots of g(x,s) in the half-interval < x ^ tt for s G (1,2). Improvement 
of estimate (25). Since g(0, s) = and 17(71", s) = — sin(7rs) > 0, we see that from existence of a root 
in the interval (0, tt) (at tt there is no root) it follows that g(x, s) in [0, tt] has a minimum at some point 
in the interval (0,7r). Derivative of g at this point is equal to zero and the value of g is nonpositive. 
Thus, existence of a root g(x, s) in the interval < x < tt for s G (1, 2) implies existence of a critical 
point in (0,7r), where the value of the function (7 is nonpositive. 

Let us find critical points (roots of the derivative by x) of the function g(x, s) for any s. We have 
g' x (x,s) = s(cosx — cos(sx)). Consequently, 

. I X = SX — 27TO1, 111 GZ 

<7a;(x, s) = <=>- cosx = COS sx < _ <^ 

I —x = sx — 27rn2, ri2 G Z 

^ xG {^l| niGZ }, xG {^| n2G Z}. (26) 
s — 1 s + 1 
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Since s G (1,2), we see that the interval (0,7r) contains only one point of form (26). Namely, 
x = Compute the value of g at this point: 

2vr \ / 2vr \ f 2tts\ ( 2tt 

s sin — sin = s sin I + 



s + 1 J \s + lj \s + lj \s + l 

2-ks \ ( 2tt \ ( 2tt 

+ sin 27T = s sm + sin 



s + lj \s + 1/ \s + l 

= (S + I)sin(-| T )>0. 

This inequality contradicts the conclusion (presented above) from the assumption of existence of a 
root g in (0, tt). Hence, g(x, s) has no root in (0, tt] for s £ (1,2) and we obtain the following inequality: 

. irm , 1 . 

xi(s) >7r for s€ (1,2) =>pi > , m€(0,-). 

1 — m 3 

So, we improve estimates (25) as follows: 

(irm \ 37rm 1 
pm, < pi(m) < , < m < -. (27) 

1 — m J 2 3 

It is important that the functions in lower and upper bound (27) have the same limits as m — > 0+ 
and as m — > | — 0. Also it is easy to study the case m = | (i.e. s = 2): 

pi G) = i 

4.4. Preliminary two-sided estimate of x\(s), s > 2. First we prove that 

7r — arcsin ( - J ^ xi (s) ^ 7r + arcsin ( - ) Vs > 2. (28) 



Existence of a root on the interval tt — arcsin Q) ^ x ^ tt + arcsin Q) immediately follows from the 
relations 

g ^7r ± arcsin ^-^ , s^j = ±1 — sin ^tts =F s arcsin ^~ 

g (^k — arcsin ^-^ , s^j ^ 0, g ^7r + arcsin ^-^ , s^j ^ 0. (29) 

Therefore x\(s) ^ tt + arcsin (i). To prove estimate (28) it remains to check that for any s > 2 
the function g(x,s) has no roots for < x < tt — arcsin Q). It is clear that ssinx > 1 when 
arcsin (i) < x < tt — arcsin (-). Hence, g(x, s) > 0. In the half-interval < x ^ arcsin (-) the 
function <?(x, s) increases (due to decrease of cost on ^ t ^ tt we have cos(sx) < cosx for < x < 
arcsin (i) =^ g' x (x,s) = s(cosx — cos(sx)) > 0). Since g(0, s) = 0, we see that the function <?(x,s) is 
positive for < x ^ arcsin Q). So, two-sided inequality (28) is proved. Now we improve it. 

Lemma 4.1. If s G (2k, 2k + 1), k G N, then the function g(x,s) has a unique root in the half-interval 
tt — arcsin Q) ^ x < tt. If s 6 (2k + 1, 2k + 2), /c G N, £/ien t/ie function g(x, s) has no roots on the 
interval tt — arcsin Q) ^ x ^ 7r and /ias a unique root in the half-interval tt < x ^ 7r + arcsin Q) . 

Corollary 4.1. The function xi(s) has the following estimates 

(tt - arcsin (±) ^ Xl (s) < tt, se(2k,2k + l), 

[tt < xi (s) ^ 7T + arcsin (±) , s G (2/c + 1, 2k + 2), G N. 
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Proof of Lemma J^.l. Assume that s G (2k, 2k + 1), where G N. Then, since g(-rr, s) = — sin(7rs) < 
and g (tt — arcsin Q) , s) ^ 0, it follows that there exists a root of g(x,s) in the half-interval tt — 
arcsin ^ x < tt. Uniqueness of the root follows from the constant sign of the derivative g' x (x,s) 
that is equivalent to absence of its roots. Let us show that for all s > 2 and x G I s = [tt — j^,tt] we 
have g' x (x, s) / 0. Notice that I s contains the interval [tt — arcsin Q) , tt], since for any s > 1 we have 
s arcsin Q) < |. For all 2k < s < 2k + 1 we see that I s does not contain any critical point (root of 
g' x ) presented in (26), since 

27r(fc-l) tt 2-rrk 2irk tt 2vr(/c + l) 

" r~ L <' TT -—, 7 > TT, 7<K-—, " ; — > TT. 

s-l 2s' s-l s + 1 2s' s + 1 

The second and the fourth inequalities are absolutely obvious. Let us check the first and the third 
inequalities. We have 

2-rr(k -1) tt 2k -2 1 , s-l , 11 

that is true, since even s > 2k. Further 

2-rrk tt 2k 1 , s + 1 ,11 
< tt - — 44> < 1 - — 44> 2k < s + 1 — 44> 2k < s + - - 



s + 1 2s s + 1 2s 2s 2 2s 

The last inequality is also true, since s > 2k and i — ^- > 0. 

Now, assume that s G (2k + 1,2k + 2), where k G N. Let us prove that g(x, s) has no roots on I s . 
We have 

g (^-rr — — , s^j = s sin — sin (ns — — ^ = s sin (t^ - ) + cos(-7rs) > 1 + cos(-7rs) > 0, 

g(ir, s) = — sin(7rs) > 0. 

Thus, existence of a root on I s would imply existence of a critical point x G I s such that g(x, s) ^ 0. 
It is easy to prove that I s contains only two critical points (26), namely 

2Trk 2ir(k + 1) 

s-l' s + 1 ' 

Let us calculate the values g (^j+\, s^j , g ^ 27r ^ 1 - 1 , s^j . We have 

2-rrk \ ( 2-rrk \ ( 2-rrk s 
, s = s sm — sin 



s-l 

/ 2-rrk \ ( , 2-rrk \ . f 2-rrk , 
s sm - sin 2-rrk -\ = (s - 1) sm > 0, 



/2vr(A; + l) \ . /2vr(A; + l)\ . /2tt(£; + 1)s 
g\ — -,s =ssm — - sin 1 



s+1 J \ s+1 J \ s+1 



= s sm 



(^)-*(^ + .)-«t±2)-(. + i)-.(*S±2)>a 

Thus we see that g(x, s) has no roots in I s . Then, since g(-rr, s) = — sin(7rs) > and <7(7r+arcsin ^) ^ 0, 
it follows that there exists a root of g(x, s) in the half-interval tt < x ^ tt + arcsin Q). Uniqueness of 
the root follows from the inequality g' x (x, s) ^ for all s > 2 and x G I s = [tt, tt + ^] . Note that I s 
contains the interval tt ^ a+(s) ^ tt + arcsin Q). For all 2k + 1 < s < 2k + 2 we see that I s does not 
contain any critical point presented in (26), since 

2-jrk 2-K(k + l) tt 2rr(k + 1) 2rr(k + 2) tt 

< TT, : > TT + —. < TT, > TT + — . 



s-l ' s-l 2s' s + 1 s + 1 2s 
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The first and the third inequalities are absolutely obvious. Let us check the second and the fourth 
inequalities. We have 

2vr(fc + l) tt 2A; + 2 i 1 n , n , s - 1 n , n s + 1 

— i — -^>7r + — & ->l + _^ 2k + 2>s-l + — — & 2k + 2>s-—- , 

s — 1 2s s — 1 2s 2s 2s 

that is true, since > and even s < 2k + 2. Further 

2tt(A; + 2) tt „, s + 1 „, „ s + 1 „ 1-s 

— i ^>7T + — 2/c + 4>s + l + 2£; + 2>s-l + 2 A; + 2 > s + . 

s + 1 2s 2s 2s 2s 

1-8 

2s 



The last inequality is also true, since s < 2k + 2 and 4^ < 0. This completes the proof of Lemma 4.1. 

□ 

Proposition 4.1. There holds the equality 

xi(s) = vr, VsGN, s^2. (31) 

Proof. Let us prove (31). Since s G N it can easily be checked that tt is a root of g(x,s). In view of 
(28), it remains to prove that g(x,s) has no roots in the half-interval tt — arcsin Q) ^ x < tt. Put 
y = tt — x. Thus we have 

g(x, s) = g(ir — y,s) = s sin(-7r — y) — sin(-7rs — ys) = s siny — 
— (— l) s sin(— ys) = s siny + (— l) s sinys ^ s siny — sinys = s). 

It was proved above that g(y, s) > for y G (0, arcsin ^] . So, s) > for 7r — arcsin Q) ^ x < tt. 
This completes the proof of Proposition 4.1. □ 

Remark 4.3. Note that inequalities (30) turn into equalities in an infinite set of values of the param- 
eter s. In fact for any k G N there exist s\^ G {2k + ^, 2k + l) and S2,fc £ (2A; + 1, 2k + |) sitc/i i/iai 

xi{si,k) = tt - arcsin (577) xi(s 2 ,k) = tt + arcsin • 

4.5. More accurate estimate of x±(s). Let p(z), 2 G M. be the distance from z to the closest 
integer number. For any s ^ 2 define a function o(s) as follows: 



a(s) 



mm 



(fp(f), arcsin (I)) for < p (§ ) < 
min (f (1 - 2p (f ))* , arcsin for ^ p (§) < I. 



Note that a(s) is nonnegative, is equal to zero for all integers s, and tends to zero as s 4 +00. 
Absolute value of derivative a'(s) is equal to ^ when s is odd number and is equal to 00 when s is 
even number. 

Proposition 4.2. x\(s) admits the following estimates: 

tt - a(s) ^ xi(s) < tt, 2k < s < 2k + 1, 

tt < xi(s) < TT + a(s), 2k + l<s<2k + 2, k G N. 

Proof. First assume that s G (2fc, 2k + 1), where fc G N. In view of (30) and (28) it remains to prove 
that 



'vr - f p (f ) <xi(s), sG(2fc,2fc + |), 
tt- fts <xi(s), s G (2/c + X, 2/s + 1), 



(32) 



where t = 1 - 2p (|) = 1 - 2 (| - fc) = 2fc + 1 - s. (33) 
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It follows from Lemma 4.1 that g(x, s) has a unique root in the interval < x < n. Since g(ir, s) < 
we see that to prove (32) we must prove that 

9(*-^p{%)> 8 ) >0 ' se(2k,2k + ^Y (34) 



- I \ n 

9 I 7T — ,s I > 0, s £ 



2ife + |,2fc + l). (35) 



We have 

* (* " 7 P (I) ' s ) = ssin (* " 7 P (I)) " sin (" s " * p (I)) • 

Note that s = 2k + 2p (§) for s £ (2k, 2k + 1). Therefore, we have 

~ 9 ( n ~ ~ s p (I) ' S ) = SSin G P (I)) " Sin H + * p GD) = SSin (I)) " Sin (I)) > °- 

(We already noticed that sin a < nsin (^) for a G (0,7r), n > 1.) Thus inequality (34) is proved. Let 
us prove (35). We have 

(irts \ ( irt^\ ( i\ 

7r , s I = s sin I 7T I — sin y-Ks — irt 'i J . 

Using (33), we get 

gi^Ti — ^— , = s sin ~~ s ™ + 1 — — 7rt 3 ) = 

^— I — sin (2irk + it — lit — nt 3^ = s sin ( ^— ■ J — sin (jrt + 7rf3^ . 



s sin 



s I \ J \ s 

Therefore, we must prove the following inequality: 



7T^ \ „ „ 1 _ 1 



sin (nt + ttU^ < ssin ( - — j for < i ^ -, s^3 



Note that vrt + Tr^ vr (§ + §) = ^ < 2 and if < a < 2, then sina <a~if + m<a-^- + j^ 



3 „ ., ., , . „,3 



a — =f|-. On the other hand, sina > a — for any a > 0. Hence, we need prove that 



i\3 fnU ttH\ ttH 2 / i\ 3 



The right-hand side of this inequality exceeds ^vr 3 t. Therefore, the required inequality will be proved 
if we verify that 

i ^ 2 2 ,r 2 A, 2 - 5 ' 

1+ 6^ < 15 7r ^ 15<7r 



Since s > 2.5, we have 
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Now consider s G (2k + 1,2k + 2), where k G N. It is sufficient to prove that g(7r + a(s),s) ^ 0. 
According to the definition of a(s) this follows from the two inequalities 
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(tt + -U , sj < 0, s£ ^2A; + 1, 2fc + 
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t = s-(2fc + l), (36) 



5(vr+%(0 ,s) <0, S £^+^2H2j. (37) 

(The inequality g (tt + arcsin Q) , s) ^ was proved above.) 

In fact, inequalities (36), (37) are in some sense "symmetric reflections" of inequalities (34), (35) 
and follows from them. To verify this we prove the following implication (with account of the following 
lemma, the proof of Proposition 4.2 is complete). □ 

Lemma 4.2. Let 5 G (0, § ] , s_, s + G R, 1 < s_ < s + , s_ + s + = 2u, where u G N. // 
5 (tt - £, a-) > 0, i/ien g (n + 1 s+) < 0. 

It is easy to show that inequalities (34), (37) and (35), (36) form exactly such pairs, as in Lemma 4.2. 
According to this we have that (34) implies (37), and (35) implies (36). 

Proof of Lemma 4.2. It is sufficient to prove that the sum g (tt — j-, s-^j +g (tt + is negative. 

We have 

s_ sin (tt — | — sin (-7rs_ — 5) + s + sin (it + — ] — sin (tts + + 5) = 

V »-) V s + J 

s- sin ~~ s + sm ^ ~~ ( s * n ^ 7rs_ — <5) + sm ( 7rs + + <$)) • (38) 

If S G (0, |] , then the function ssin (|) increases for 1 < s < +oo and sin _ s + sin < 

0. Since sin ^ 7r ( s -+ s +) ^ _ s i n (7^), 1/ g N, we have (sin (7rs_ — 5) + sin (7rs + + 5)) = 0. This completes 
the proof of Lemma 4.2. □ 

Figure 2 shows the plot of x\(m) and obtained two-sided estimates (dashed lines). Since the function 
decreases in the interval < t < tt, we have arcsin Q)<^fors>2. Therefore 



tto fl\ „ 1 

— ^ arcsin I - , if p ^ 



3' 



-{/l-2p^ arcsin ( - j , if ^/l - 2p ^ i ^ p ^ H 
s V s / 3 27 



So if p G [|, ^|], then arcsin is the best estimate. Therefore in Theorem 4.1 the minimum is a 
well-posed operation. Returning to the original variables p = x^-, m = we obtain the statement 
of Theorem 4.1. 

4.6. Study of p\(m) for m > 1. We use the change of variables p = ^, m = ^ G (0, 1) to study 
Pi(m) in the case m > 1. We have g\(p,fh) = — ^gi(p,m). Therefore, we have g\(p,rh) = iff 
gi(p, m) = 0. In such a way we get the following functional equation: 

Pi(m) = mpi( — ). (39) 
m 

Hence, the properties a) — f) (see Theorem 4.1) of p\(m) for m > 1 follow from the similar properties 
of pi(m) for m G (0, 1). 



12 



Fig. 2. Plot of x±(s) with two-sided estimates (dashed lines) 



4.7. Differentiability of pi(m). In this subsection we prove that the function p\ (m) is a continuous 
and decreasing function for all m > 1 and pi(m) is continuously differentiable at all points of the 
interval m E (1, +oo) except for the set {m = ^-^\k E N} where its derivative is equal to — oo. From 
functional equation (39) it immediately follows that the function p\ (m) is a continuous and increasing 
function for all m E (0, 1) andpi(m) is continuously differentiable at all points of the interval m E (0, 1) 
except for the set {m = ^i\k E N} where its derivative is equal to +oo. 
Let m > 1. From (31) and (39) it follows that 

k + 1 

pi( — - — ) = ir(k + 1) for any k E N. 

k 

Consider the partial derivatives 

dgi{p,m) m 2 -l . . p 

sin p sin — , (40) 



dp mm 



dgi(p, m) _ mpcos(p) cos (£) + (psin(p) - m 2 cos(p)) sin v? 



, > ( 41 ) 

We have ^(tt(A; + 1), = and §^(tt(A; + 1), = -ttA; < 0. Hence the plot of p x (m) has a 

fc+i i 



vertical tangent line. We claim that pi(m) is smooth for m > 1, m ^ {^r- |& E N}. Indeed, this follows 
from the implicit function theorem and the inequality ^-| p=pi ( m ) ^ 0. Let us prove the last inequality. 

The equation = has two series of positive roots, namely p l {k\) = 7rfci and p 2 (k2) = irk2m where 
k\,ki E N. At these points we have 

gi(p 1 (ki), m) = (— l) fel+1 msin(— !— ), gi(p 2 (k2),m) = (— l) fc2 msin(/c 2 m7r). 

m 
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It can easily be checked that gi(p l (ki),m) / for m G (f+fj^rO =: ^(&)> Vfc G N. This implies 
that ^-| p=pi ( m ) 7^ 0. From the implicit function theorem it follows that pi(m) is a continuously 
differentiable function on I(k) and 

dgi (p, m) dg\ (p,m) p + m cot p{—m + p cot ^) 



PiH = ^ — / 



p=Pi(m) — ™f™2 _ "H lp=pi(m)- 



5m dp ( ' m(m 2 — 1 

Now we claim that p\{m) monotonically decreases for m G (1,2), to ^ {^|r-|& G N}. Indeed, this 
follows from the inequality p'^m) < 0. Let us prove this inequality. Since gi(p, m)| p=pi ( m ) = 0, 
we get TOcotp| p=pi(m) = cot ^| p=pi(m) . Then sign(pi(m)) = -sign(/(pi(m), m)) where f{p,m) = 
p — m cot — + pcot 2 — . Consider the function / as a polynomial of second degree with respect to 
cot — . Since to < 2 and pi(m) > p > 1 (see Theorem 4.1 e)), we see that the discriminant of the 
polynomial is negative, indeed D(/) = to 2 — \p 2 < 0. Thus f{p\{m),m) > and Pi(m) < 0. 

Now we claim that the function pi(m) is continuous at the points m = for any fceN. Indeed, 
this means existence of limits of pi(m) at the considered points and lim pi(m) = pi(^r-) = 

ir(k + 1). Since pi(m) is monotonic and bounded for to / ^M- it follows that there exists finite 
limits p±(k) = lim p\{m). Inequalities p+{k) < 7r(/c + 1), p~(k) < ir(k + 1) contradict to 



k 

'k+l 



*±1±0 



Pi{—%-) = Tt(k + 1) because if they are satisfied, then there exists a positive root p < pi(m) of the 
function gi(p, m), but by definition p\(m) is the minimal positive root. In such a way we proved that 
p+(k) ^ ir(k + 1), p-(k) ^ n{k + 1). Further, the inequalities p+(k) > ir(k + 1), p-(k) > ir(k + 1) 
contradict to continuity of the curve {(p, m)\gi(p, to) = 0} in a neig hborhood of (ir(k + 1), but 
from the implicit function theorem this curve is continuous because ^^{^{k + 1), = — m ^TT 7^ 0- 
Thus we proved that p+(k) = n{k + 1), P-(k) = ir(k + 1). So p\{m) is a continuous function for any 
m G (1,2]. 

Now we show that Pi(^-r^) = — oo for any k G N. For any to* = we know the explicit value 
P* = Pi{rn*) = ir(k + 1). We have d91 %^ n - = —Ttk and dgi(p . m ) _ q This implies existence of the 
limit lim p'Am) = — lim jp-/^f-^-(pi('rn),m) = — oo. Therefore there exists p\ (to* ) = — oo. 

Thus we proved that pi(m) is continuous for any to G (1,2] and p\ G C 1 (il) where := {to G 
(1, 2)|to / ^±i, fe G N}. Now let to > 2. We have ^-| P=Pl ( m ) + 0. We claim that p[(m) < 0. Indeed, 
this means that f(p,m) > for p = pi(m). In fact, if D(/) < 0, then f(p,m) > 0; if D(/) ^ 0, then 
m ^ 2p and we get cot ^ > <?2 where q2 is the greatest root of the function / as a polynomial of second 
degree with respect to cot — . Thus, we proved that f(pi(m),m) > 0. It follows that p'i(m) < 0. So 
pi(m) is a continuously differentiable decreasing function for to G (2,+oo) 

Finally we show that the function p i— > g\ (p, to) changes sign at the point p = p\ (to) . If m £ { \ k G 
N} this follows from the inequality ^-| p=Pl ( m ) / 0. If to = to* = we have ^-| p=pi( - m ) = 0, 

^jr\p= Pl ( m ) = 0, ^| p=pi ( m ) = - 2 ^+iy2 < 0. It follows that gi(p,m*) changes sign at the point 
P* = Pi{m*). 

5. Study of the function P2{m) 

In this section we prove a two-sided estimate of a minimal positive root of the equation 52 (p, m) = 0, 
where 

92{p, m ) = mp cos (—) sinp — (pcosp + (to 2 — 1) ship) sin ( — J . (42) 

Our aim is to find or estimate as accurately as possible for any m G (0, 1) U (1, +00) the value of P2(rn) 
defined in (17). The main results on this problem are summarized in Theorem 5.1. 
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Theorem 5.1. For any m G (0, 1) U (1, +oo) the function giip-, m) has a minimal positive root P2(m) 
satisfying the following properties: 

a) The function p2(m) is continuous for m G (0, 1), coincides with at the points {m = ^j\k G 
N} and at points m = m\ k , where — 15+ 40fc + 32fc2+8fc 3 < 171 1 k < 3+2fc > ^ continuous for 
m G (l,+oo), coincides with at the points {m = ^-\k G N} and at points m = m* 2k , 
where f±|| < m\ k < f±|| + 1+gA , 2 1+A ,^ • 

b) T/ie function P2{m) is continuously differentiable at all points of the set (1, ^) U (2,+oo) and 
the following sets: 



m G 



k l + 2k 



1 + k ' 3 + 2k 15 + 40k + 32/c 2 + 8/c 3 



U 



U 



/c + 2 3 + 2/c 



fc + 1' l + 2k 



u 



3 + 2A; 



+ 



l+2fc 1+jfc' 
3 + 2/c' 2 + /c 
2 fc + 1 



U 



1 + 2A; 1 + 8fc(l + k 2 ) ' fc 



/or aZZ k G N. 0w£ o/ i/iese sets i/tere exist values m = m\j~, g+fjf — is+40fc+32fc 2 +8fc 3 ^ m i fc < 
m hk < 3+H> anrf m = »™2,fc, f+H < ™2,fc < m* 2 k < f±|| + 1+8fc ( 2 1+ ^) , ^/iere i/ie derivative 
p' 2 (m) is equal to +oo and — oo respectively. 
c) For m G (0, 1) the function P2(m) has the following lower and upper bounds: 



5.7m < P2{fn) < 2nm for m < ^, 

(S)<l*(m)<S + s£lA™>3- 



7rm _ m 
1— m 1— m 



arcsm 



(43) 



For m > 1 £/ie function p 2 (rn) has the following lower and upper bounds: 



-™l _ ™. arcsin ( a=i ] < p t m ) < 4 

m— 1 m— 1 \ m+1 / ' m— 1 



3— m 



for m ^ 2, 



[5.7 < P2{m) < 2tt for m > 2. 
d) For m G (5, 1) pZoi 0/ p = p 2 (m) "wraps" the hyperbole p = as follows: 

P2(m) <^ for me ±±f - 15+40fc+ 2 32 ^ +8fca 



(44) 



P2{m) > f=^i f or m G 



l+2k l+k 



k G N. 



(45) 



(46) 



3+2fc> 2+fc / 

For m G (1,2) the plot of p = pi(m) "wraps" the hyperbole p = as follows: 
( P2(m)>^/ rmG(|±f,f±f" , 

' P2(m) < ^ /or m G [f±§ + 1+8fc(1+fc . } , *p) , VA; G N. 

e) At the point p = P2^jn) the function p 1— > g2(p, m) changes its sign. 

Figure 3 shows the plot of P2 (m) with two-sided estimates (dashed lines) . 

Remark 5.1. From (43) and (44) it follows that p2(m) tends to +00 as m — > 1 ±0. 

5.1. Two-sided estimate of p 2 {m) for < m < \. Let us prove that for m G (0, |) the equation 
g 2 (p,m) = (see (42)) has a root in the interval < p < 2irm. First write the equation g 2 (p,m) = 
in the equivalent form 



mpcos ( — ) ship = (pcosp + (m 2 — 1) ship) sin ( — ) 
\mJ ' \m) 



(47) 
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Fig. 3. Plot of the function p2(m) with two-sided estimates 



Since < m < ^, we see that in the interval < p < 2irm the product sin p sin (^-) vanishes only at 
p = irm but it can easily be checked that this point is not a root of equation (47). Therefore, dividing 
both sides of (47) by the product sinpsin(^-), we get the equivalent equation 

mpcoti— ) = pcotp + m 2 — 1 44> mpcoti— ) — m 2 = pcotp — 1 44> 
m m 

mpcot(^) — m 2 pcotp— 1 ^ cot(^) — 1 _ pcotp — 1 

p 2 p 2 (£_j2 p 2 

Let G(x) = 1 ~ x ^ tx . We have 

V 

g 2 (p, m) = 44> G(-) = G(p) for < p < 2nm. 

m 

Put x = Rewrite the last equation as 

G(x) = G{mx), < x < 2vr. (48) 

Since the function G(x) decreases over the intervals < x < it and ir < x < 2ir, we see that equation 
(48) has no roots in the interval < x < ir. Moreover, since lirriaj^o G(x) = |, it has no roots in the 
interval it < x < p%, where p2 is the root of the equation G(x) = | (p2 = 5.7634. . . ). On the other 
hand, from the same reasons related to monotonicity of G it follows that for any x £ (p2,2ir) there 
exists a unique value m = m{x) G (0, |) such that there holds equality (48). 
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We claim that this dependence m(x) is differentiable and has a positive derivative. Indeed, we must 
prove that 

, G'(x) - m(x)G' (m(x)x) 

m (x) = . . r . > 0. 49 

xG'(m(x)x) 

Inequality (49) follows from the following inequality: 

G'(x) - mG'(mx) > for G(mx) - G{x) = 0. 

Let us prove it. We have 

1 cotxV l + cot 2 x 1 G(x) 



G'{x) 



rp2i ryt f rp r f 3 ryt 

tXi if,' / tXi tXi if,' 



mG'(mx) 



m 



l + cot 2 (mx) 1 G(mx)\ l + cot 2 (mx) 1 G(mx) 



mx (ira) 3 mx ) x m 2 x 3 x ' 

where G'(mx) is the derivative G'{x) taken at mx. Since G ^ x ) = ^iM we must prove the following 
inequality: 

cot 2 (mx) cot 2 (x) 1 1 i, \ i, \ 1 1 
- < — =— r — — ^ <^ cot (mx) - cot (x) < — =— ~ — 

x x m z x a x 6 m z x z x z 

Since coti < \ for any t G (0,7r), we have cot 2 (mx) < (mx)~ 2 . Further, since for any x G (p2,2tt) we 
have cot x < — 1, it follows that cot 2 x — x~ 2 > 1 — p^ 2 > 0. Thus, we proved the required inequality. 

Now, since the dependence m(x) is differentiable and has a positive derivative we can use the inverse 
function theorem and conclude that for any m G (0, ^) there exists the differentiable and increasing 
function x(m) G (p2, 2tt) such that x(m) is a root of equation (48). Thus, we see that P2(m) = mx(m) 
is a differentiable increasing function for m G (0, ^), which admits the following estimate: 

Pim < P2(m) < 2-jrm, < m < —. (50) 

In addition, we claim that for any m G (0, ^) the function P2(m) exceeds both po(m) = an d 
Pi(m) (see Theorem 4.1). Indeed, for < m ^ | this follows from the inequalities ^ P\{m) ^ 
and P2(m) > p2m > ^!p. For | < m < \ we have p\[m) ^ and it remains to prove that 

irm . , 1 1 , ri . 

< P2{fn), - < m < -. (51) 



1-m ' 3 2 

In fact, it can easily be checked that for any m / 1 there holds the equality cot( 33^) = cot( and 
for any m G (|, |) there holds the inequality cot(j^) < 0. By definition, put A = — cot(j^) = 
— c °t(l?^) an d Gi(x) = G(x) — G(mx). The function Gi(x) increases over the interval -it < x < 2tt. 
It follows from definition of p 2 {m) that G^^M) = 0. Since G\(x) increases, we have 

(51) * Gi(-^-)<0 ^ G (-^)<G(-^-) ^ 
1 — m I — m 1 — m 

1-,-*- cot (,-*-) l-^-cot(^) 

/ it \2 ^ / 7rm \2 

7T.4 1 + tt^ 1 fi 1 



1 — m m 2 1 — m m 2 1 — mm 

The last inequality is obvious since 1 < < i?4im- Thus, inequality (51) is proved. 
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5.2. Reduction of (72 (p,m) = to a simpler form. In this subsection we reduce the problem of 
finding (estimation) of P2(jn) for m 6 [5, 1) to finding (estimation) of the minimal positive root X2(s) 
of a simpler equation #2(x, s) = for s ^ 3, where 52(2^, s) = 4s(cosx — cos(sx)) — x(s 2 — l)(ssinx + 
sin(sx)). 

It can easily be checked that 



. p m + 1 p 

92\P, m ) = — ym — 1) sinpsin — ) + p( fsmpcosf — ) — 

m 2 m 

p — \ p p 

— cospsin( — )) H (sinpcos( — ) + cospsin( — ))). 

m 2 mm 

Transform the first summand: 

(m 2 - l) sinpsin (£) = (cospcos (£) + sinpsin (^) - cospcos (£) + sinpsin (£)) = 

= ^^^l^ («* (P (1 " £)) " COS (p (1 + i))) = 

= 4(1 + SS)" 2 *«S» (cos - cos (S^)) ■ 

Further, transform the second summand: 

V sin (p - £) + V sin (p + £)) = f (m - 1) sin (p - £) + sin (^^) ) = 

= K 1+ ^)" (- 1 + (^) 2 ) (^) sin (^) + sin (^^)) • 

Let us introduce a parameter s = y±^- and a new variable x = P^p- Now express g2(p,m): 

§2{p, to) = 2(i+ s ) i 9 2 ( x > s )' wn ere £b(x, s) = 4s (cos x — cos (sx)) — x (s 2 — l) (s sinx + sin (sx)) . 

Thus, we reduce the problem of finding (estimation) of P2(m) for m G (5,1) to finding (estimation) 
of the following function: 

x 2( s ) = min{x > 0\g2(x, s) = 0} for s ^ 3. 

In Subsections 5.3 - 5.4 we obtain a two-sided estimate of X2(s). Returning to the original variables 
p = x^-, m = J=j we obtain the statement of Theorem 5.1. 

5.3. Two-sided estimate of X2(s) for s ^ 3. In this subsection we prove that X2(s) admits the 
following lower and upper bounds: 

7r — arcsin - < X2(s) < ir -\ — for s ^ 3, (52) 

s s — 2 

and X2(s) is a unique root of 52 (x, s) on the interval (it — arcsin 7, vr + 5^2). 
Proposition 5.1. There holds the inequality 

3-7T 

x 2 (s) < — . (53) 

Proof. Indeed, (53) follows from the fact that 32 s) has zero derivatives up to the fifth order at the 
point x = and the sixth derivative is negative. In fact, - g^e'^ U=o = — 2s(s 2 — l) 3 < for all s > 1. 
It means that the function g2(x, s) is negative in some right half-neighborhood of the point x = 0. 
But at the point x = ^ this function is positive, since 

a,(£,.) = -4»cos (*£) - I - 1) (sin (^) -.) > -4- * ( 5 * -!)(!-.) = 

= y f s 3 - s 2 - s (l + A) +l) > s 2 - s - 2 > for any s ^ 3. 
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Hence, we see that the continuous function c/2(x,s) changes its sign over the interval (0, This 
implies that g2(x,s) has a root in the interval (0, -^). Thus, inequality (53) is proved. □ 

Now we rewrite gi (x, s) in the following form: 

c/2(x, s) = As (cos x — cos (sx)) + xH(x, s), 

where H(x, s) = — (s 2 — l)h(x, s) and h(x, s) = ssinx + sin (sx). 
Let J(x, s) = ffjf^y- It can easily be shown that 

dJ(x,s) g(x,s) 2 

= — -77 ^ 0, where g(x, s) was denned in 4.2. 

ox h(x,s)^ 

Thus we can see that if g(x, s) and h(x, s) do not vanish at the same time (this case is considered 
separately in Remark 5.2), then J(x,s) increases at all points, where h(x,s) / 0, and has vertical 
asymptotes at the points x(s) = {x\h(x,s) = 0}. The function J(x,s) has the following properties: 
1) lim x ^ J(x,s) = 0, 2)lim^ 5 , fc ( s) _ J(x,s) = +oo, 3)lim Mit(s)+0 J(x,s) = -oo, where x k (s) is the 
fc-th positive root of h(x,s). Thus, we conclude that the function J(x,s) (and hence the function 
92(x,s)) has no roots in the interval (0, xi(s)) but it has a unique root in the interval (xi(s), X2(s)). 
Hence, we arrive to the problem of finding (estimation) of x\(s) and X2(s). We study this problem 
in Proposition 5.2 (see below). Combining the estimations X2(s) G (0, ^-) and x\(s) < X2(s) < X2(s) 
with Proposition 5.2 we get the following two-sided estimation of X2(s): 

7r — arcsin - < X2(s) < (54) 

Remark 5.2. Consider the case when g2(x, s) and H(x, s) have a common root for < x < This 
holds if and only if h(x, s) = and cosx — cossx = 0. In Subsection 4-3 (see (26)) we found the roots 
of the second equation, namely x G {|^j|n G N}. Let us find the values of s, for which these roots 

are the roots ofh(x,s) simultaneously, in other words ssin|g- + sin (^Sj^J = 0. The last equation is 
equivalent to (s =F l)sin||j = <^ -Jg- = k, k G N. Thus, s = ^ -f 1. So, we got that g2{x,s) and 
H(x,s) vanish simultaneously if and only if x = irk and s = -|- =p 1, where k, n G N. It is clear that 
only for k = 1 and s = 2n + 1 this solution satisfies the limitations x G (0, -jj-) and s ^ 3. Thus, we 
get 

g 2 (x, s) = H(x, s) = O x = 7T, s = 2n + 1, n G N. 

Since J(0, s) = 0, J(x, s) increases over the interval x G (0, it), and h(x, s) > for s = 2n + 1 in this 
interval (see Proposition 5.2), we see that J(x,s) > in this interval. This implies that for s = 2n + l 
the function g2(x, s) has no roots for x G (0, tt) and x = n is the minimal positive root. Thus, we get 

x 2 (s) = tt for s = 2n + 1, n G N. (55) 

In the following proposition we estimate the first and the second positive roots of the function 
h(x, s) = ssinx + sin (sx). 

Proposition 5.2. The functions 

x\(s) = minja; > 0\h(x,s) = 0} 

and 

X2(s) = min{x > x\(s)\h(x, s) = 0} 
have the following two-sided estimates for s ^ 3: 

~ / x r ■ 1 , 11 rl 

xi(s) G tt — arcsin -; tt + arcsin - =L, 
s s 
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1 1 

2tt — arcsin - ; 2tt + arcsin - 

s s 



I 2 



x 2 (s) G 

This estimate for xi(s) can be improved as follows: 
x\(s) = it, for s£N, 

xi(s) G [it- arcsin ±; ir) , s G (2k + 1,2k + 2), 
xi(s) G (vr;7r + arcsin ±] , s G (2fc + 2, 2k + 3), fc G N. 

Proof. First, we show that there exist roots of h(x, s) in ij and if 

f 1 \ ( 1 

n 7r =F arcsin -, s ) = ±1 + sin [it =f s arcsin - 

V s J V s 

h (^k — arcsin -, s^j ^ 0, h ^ir + arcsin -, s^j ^ 0, (56) 

h ^2ir =F arcsin -, s^j = =Fl + sin ^2irs =F s arcsin - 

h ^2ir — arcsin -, s^j ^ 0, h ^2tt + arcsin -, ^ 0. (57) 

Further, if arcsin (^) < x < 7r — arcsin i, then s) > because in this interval we have ssinx > 1. 
Now we claim that h(x,s) > for < x ^ arcsin^. Indeed, this follows since h(x,s) is increasing 
over this half-interval and from equality h(0, s) = 0. The last equality can easily be checked by direct 
calculation. Let us prove that h(x, s) increases over the half-interval x G (0, arcsin ^] . To prove this 
we consider the derivative h' s (x, s) = s(cos x + cos sx) and find all positive critical points 

. I x = sx + 7r — 2-irni , n\ G N 

h x (x, s) = <^ cosx = — cossx 44> < <^ 

I —x = sx + 7T — 27rn 2 , n 2 G N 

x G {^rvr|n G N}. (58) 

Since for any s ) 3 we have minj^gp7rj = > ^ > arcsin -j:, we see that the half-interval 

x G (0, arcsin i] does not contain any point of form (58). Moreover, h' x (0, s) = 2s > and we proved 
that h(x,s) is increasing. 

So we proved that h(x, s) > for < x < tt — arcsin i. Now we claim that there exists a unique 
root, denoted by xi(s), of the function h(x, s) in the interval tt — arcsin ^ ^ x ^ tt + arcsin ^. To prove 
this we divide the set s ^ 3 into the three subsets s = k+2, s G (2/c + l, 2k + 2), and s G (2k + 2, 2k + 3), 
where k G N. 

(1) Let s = k + 2, k G N. We have /i(-7r, s) = 0. Uniqueness of the root x\(s) = tt for tt — arcsin ^ ^ 
x ^ it + arcsin ^ follows from monotonicity of /i(x, s) w.r.t. x in the intervals it — ^ < x < it 
and 7r < x < 7r + ^. Indeed, if the first interval contained a critical point (58) then there would 

exist n G N such that (fc + 2) ± 1 — (| ± 2( ^ +i) ) < 2n - 1 < (k + 2) ± 1, but this contradicts 
to the following obvious inequality: 

WeN:0< G ± 2(Fr2)) <1 - < 59 > 

If the interval (it; it + ^) contained a critical point (58) then there would exist n G N such that 
(fc + 2) ± 1 < 2n - 1 < (A; + 2) ± 1 + (± ± 2(FP2)) , but this contradicts (59). 



20 



(2) Let s 6 (2k + l,2k + 2), fceN. We have /i(7r, s) = sin7rs < 0. Combining this with (56) we see 
that the function h(x, s) has a root in the half-interval tt — arcsin (^) ^ x < tt. Uniqueness of the 
root xi(s) follows from the constant sign of the derivative h' x (x, s) that is equivalent to absence 
of its roots. Indeed, let us show that the interval [tt — ^, tt) D [V — arcsin ^, 7r) does not contain 
any critical point (58). Assume the converse, then there exists n 6 N such that 1— ^ < < 1. 

Since 2/c + l < s < 2k + 2 we have (2& + 1 + 1) (l - 2(2fc 1 +1) ) < 2n-l < 2& + 2 + 1. Further, since 

< 2(2fc+i) = \ ^ 2(2fc+i) < 2n_1 G Nand2/c + l±l G Nwehave2n-1 = 2fc+l±l. Now we 
see that no number n satisfies the last equality, since 2n — 1 is odd and 2k + 1 ± 1 is even. This 
contradiction proves that the function h(x, s) has a unique root for tt— arcsin ^ ^ x < tt. Further, 
we show that h(-rr, s) < for tt ^ x ^ tt + arcsin (i) . It is sufficient to show that h(x, s) < at 
the end points of the considered segment and at all critical points inside it. At the left end we 
have H(tt, s) < 0. At the right end we have h (tt + arcsin \,s) = — 1 + sin (its + s arcsin < 0, 
since from (2k + 1)tt < sit + s arcsin ^ < (2k + 2)ir + f it follows that sin (sir + s arcsin ^) < 1. 
Now find critical points (58) lying on the interval (tt, -it + ^) D (it, tt + arcsin i) . It is easy to 
check that the inequalities 



(2k + 1 ± 1) < 2n - 1< (2k + 2 ± 1) j^l + 2 ( 2fc + ^ J . ^, e N 

are satisfied only for the value n = k + 2 if we select the sign " plus" , and for n = k + 1 if we 
select the sign "minus". Hence the considered interval contains only two critical points, namely 
and §±±vr. We can estimate the value of /i at these points as follows: 

, /2/e + 3 \ /2fc + 3\ /2fc + 3 
h tt, s = s sin + sin tts 



s + 1 ) \ s + 1 ) \ s + l 

f2k + 3 \ . „, 2A; + 3 

= s sin 7r + sin (2k + 3) 7r 7r 

\s+l J V 5+ 1 

2fc + 3 2A; + 3 tt 

= (s + 1) sin tt < 0, since tt < < tt H < 2tt, 

s + 1 s + l 2s 

2k + 1 \ . (2k + \ \ . (2k + \ 

tt, s = s sin tt + sin 7rs 



1/ V s - 1 / V s - 1 

. /2k + l \ . _ 2/c + l 
s sin I — 7T + sin (2k + 1) tt H — 7r 



n s-1 ) V a-1 

2/c + l 2/c + l TT 

= (s — 1) sin 7r < 0, since 7r < < tt H < 2-7T. 

y J s - 1 s-1 2s 

Thus we proved that h(x,s) <0for7r^a;^7r + arcsin i, therefore h has no roots in this 
interval. 

(3) Let s G (2k + 2,2k + 3), k € N. We claim that > if vr - arcsin (i) ^ a; < vr. It is 

sufficient to show that h(x, s) > at the end points of the considered segment and at all critical 
points inside it. At the left end we have h(rr — arcsin \ ,s) = 1 + sin (tts — s arcsin i) > 0, since 
from (2k + 2)tt — ^ < tts — s arcsin ^ < (2k + 3)tt it follows that sin (tts — s arcsin ^) > — 1. At 
the right end we have Ii(tt, s) = sin7rs > 0. It is easy to check that the inequalities 

(2k + 2+ 1) (l - ^ ^ < 2n - 1< 2k + 3 + 1, n, k G N 

are satisfied only for the value n = k + 2 if we select the sign " plus" , and for n = k + 1 if we 
select the sign " minus" . Hence the interval [tt — ^ , tt) D [tt — arcsin ^ ; tt) contains only two 



21 



critical points (58), namely ?r and ^t^tt. We can estimate the value of h at these points 
as follows: 

, /2k + 3 \ 2k + 3 -it 2k + 3 

h tt, s = is + 1 ) sin tt > 0, since < tt < tt, 

\ s + 1 ' J v ; s + 1 2s s + 1 ' 

/2fe + l \ . . 2fc + l vr 2/c + l 

/i tt, s = (s — 1) sin 7r > 0, since < tt < tt < tt. 

\ s - 1 ' J v ; s-1 2s s-1 

Thus we proved that h(x, s) > for it ^ x ^ 7r + arcsin ^, therefore /i(x, s) has no roots in this 
interval. On the other hand in view of (56) we see that h(x, s) has a root in the half-interval 
tt < x ^ tt + arcsin (|) . 

Uniqueness of the root x\(s) follows from the constant sign of the derivative h' x (x, s) that 
is equivalent to absence of its roots. Indeed, let us show that the interval (k;tt + t^) D 
(it, it + arcsin ^] does not contain any critical point (58). Assume the converse, then there 
exists n G N such that 1 < < 1 + ^. Since 2k + 2 < s < 2fc + 3 we have 

2A; + 2 ± 1< 2n - 1< (2k + 3 ± 1) ( 1 + — — ^ — — ) , 2n - 1 G N , k G N. 



2(2/c + 2) 

Thus we have 2n — 1 = 2/c + 3 ± 1. But the last equality is satisfied for no numbers n and 
k. Therefore such n does not exist. This contradiction proves that the function h(x, s) has a 
unique root for it < x ^ it + arcsin - . 

Thus we proved that for any s ^ 3 the function h(x, s) has a unique root on the interval tt — arcsin - ^ 
x ^ it + arcsin i. To conclude the proof, it remains to check that h(x, s) has no roots in the interval 
-it + arcsin i < x < 2ir — arcsin ^. Indeed, we see that h(x,s) < since in this interval we have 
ssinx < — 1. □ 

Now we improve the upper bound in estimation (54). 

Proposition 5.3. There holds the inequality 

x 2 (s) < TT + — — . (60) 

S 

Proof. Indeed, (60) follows from the inequality g 2 (tt + > s ) > that we prove now. Denote by x 
the value tt + • We have 

c/2 ^tt -\ -, s^j = 4s (cos x — cos (sx)) — x (s 2 — l) ^s sin ^tt H + sin (xs)^j > 



> —4s — 4s cos (sx) — x (s 2 — l) sin (sx) + xs (s 2 — l) sin ^ — 



(61) 



Combining this with the following inequalities: 



4s cos (sx) + x (s 2 — l) sin (sx) ^ \/ 16s 2 + x 2 (s 2 — l) 2 , 



ssin|— , ) > a | = i + _____ 



we obtain 



fc^-H^,.) >*(•»-!) (l + 7 l i - s ^)-4.- > /« ta » + ^- 1 )». 
Consequently it remains to prove that 

4s + Vl6* 2 + * 2 (s 2 -l) 2 <x(s 2 -l) (l + ^-^_L_). 



(62) 



22 



Dividing both sides of (62) by x (s 2 — l) , we get the equivalent inequality 



4s L 16fi2 2 s 



x(s 2 -l) y x2( s 2-l) 2 s-2 6(s-2) 3 ' 

Now use the fact that y/l + 2t < 1 + 1 (Vi > 0). This allows us to replace the inequality (63) by the 
following stronger inequality: 

4s 8s 2 s 2 . 

+ ,,, .,2 + TTT^W < ~ o • (64) 



x(s 2 -l) x 2 (s 2 -l) 2 6(a-2)3 a- 2 
Multiplying both sides of (64) by ^J 2 -, we obtain 

2s (s- 2) 4s 2 (s- 2) s 

X (S 2 - 1) + X 2 (s 2 _ 1} 2 + 12(s-2) 2 < • 

We strengthen this inequality by replacing the left-hand side of the product s (s — 2) with the greater 
value (s — 1) . Thus we must prove that 

2 (s - 1) 4s s 

^TT) + 1^TW + W^W <1 - ° 3 ' (65) 

If 3 ^ s ^ 4, then |=i ^ |, — s ^ 2 ^ 3, ^ jq, and the left-hand side of (65) is not greater than 

6 12 3 2 1 3 _ 2 1 _ 11 
5x + 16x 2 + 12 < 5 + 12 + 12 _ 5 + 3 _ 15' 
If s > 4, then - 2 < 1, s 2 < and the left-hand side of (65) is less than 

2 16 1 2 16 1 _ 3 16 
x + 2^ + 12 < 3 + 2^5 + 12-4 + 225 < °- 83 - 
Thus we proved inequality (65). Hence inequality (60) is proved. □ 

5.4. More accurate estimate of x 2 (s), s > 3. In this subsection we improve the two-sided 
estimate obtained in the previous subsection tt — arcsin Q) < x 2 (s) < ir + j^. 

Proposition 5.4. There hold the inequalities 

x 2 (s) >vr/orsG [2/c + 2,2/c + 3), k G N, (66) 

x 2 (s) <vr for sG (2k + 1,2k + 2-——], k G N. (67) 

2/e + 2 

Proof. To prove (66) it is sufficient to check that 52(71", s) < 0. First let us calculate the values of the 
function 

<? 2 (7r, s) = —4s (1 + cos (its)) — 7r (s 2 — l) sin (its) 

at the end points of the intervals s G [2k + 2,2k + 3) for all k G N 

5 2 (tt,2A; + 2) = -8s < 0, ^(tt, 2k + 3) = 0. 

Further, using the inequality 

cos(vrs)>-l, sin(vrs)>0, Vs G (2k + 2, 2/c + 3), 

we get that <? 2 (7r, s) < inside the considered intervals and thus (66) is proved. 

To prove (67) it is sufficient to check that 52(71", s) > 0. Let d — s — (2k + 1) G (0, 1 — 2 fc+2 

). We 

must prove that for all k G N there holds the inequality 

g 2 (iT,2k + 1 + d) = 4s(cos(vrcf) - 1) + (s 2 - l)vrsin(7rd) > 0. (68) 
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Let t = ird £ (0, tt — 2^2)- Since | < 7r — < 7r ' we nave 



sin t > t sin(7r 



7T 



2k + 2' 



t sin( 



7T 



2/c + 2 



)• 



Using the inequality Va £ (0, 5) : sin a > |a, we get 



t sin( 



7T 



2/c + 2' 



> i- 



7T TT 



> 



> 



2 2/c + 2 2/c + 2 s + 1 



Therefore we have sin t > . Combining this with Vt £ R : cos t > 1 — ^ , we obtain 



<7 2 (tt, 2A; + 1 + d) > -2st 2 + vr(s 2 - 1] 



s + 1 



t(-2st + tt(s - 1)). 



Note that t > and thus we must prove that — 2st + tt(s — 1) > 0. Since t < 1 and s > 3 we have 

-2a* + tt(s - 1) > (tt - 2)s - tt > 1.14s - vr > 3.42 - tt > 0. 
Thus (67) is proved. □ 

5.5. Differentiability of P2(m), < m < 1. Differentiability and monotonic increase of the 
function P2(m) for m £ (0, 5) is proved in Subsection 5.1. Let us prove differentiability of p 2 (m) 
for 



m £ 



/c 1 + 2k 



U 



1 + 2/c 1 + jfe' 
3 + 2fc' 2 + fc 



=: M. 



1 + fc' 3 + 2A; 15 + 40& + 32fc 2 + 8P 
Since p 2 (m) = 1^^2(13^) (see Subsection 5.2), this follows from differentiability of x 2 (s) for 

1 



s £ 



2k + I, 2k + 2- 



2k + 2 



U [2fc + 2, 2fc + 3] =: S. 



Let S = S\ {2k + 1}, £ N. Recall that x 2 (s) is the minimal positive root of <? 2 (x, s) = 0. Also, x 2 (s) 
is the minimal positive root of J(x,s) = ff|f^y (see Subsection 5.3) if s £ S. By estimate (54), we 
study differentiability of x 2 (s) in the interval x £ (tt — arcsin Q) , tt + j^) =: X. Using the implicit 

Aw(x,s)(l-s 2 )- 2 I 

g 2 (x,s) \x= ''21 

smooth function without singularities 



function theorem, we get x' 2 (s) = - 9 J ^' 8) / 8 J ^ \x=x 2 (s) 



2 t s \, where w(x,s) is a 



w(x, s) = x sin(x) sin(sa;)s 4 + (— cos(x) cos(sx)x + x + 2(cos(sx) — cos(x)) sin(x))s 3 + 
+ ((cos(sx) — cos(x)) sin(sx) — x sin(x) sin(sx))s 2 + (x cos(x) cos(sx) — x)s + (cos(sx) — cos(x)) sin(sx), 

and g(x,s) = sin(sx) — ssin(x) (we studied this function in Subsection 5.2). In Lemma 4.1 we 
proved that the equation g(x,s) = has a unique root, denoted by xi(s), in the interval x £ 
[tt — arcsin(i), tt + arcsin(i)] . It was also proved that g(x, s) 7^ in the interval x £ (7r+arcsin(i), 2tt— 
arcsin(i)). Hence xi(s) is a unique root of g(x, s) for x £ X. Thus the function x 2 (s) is differentiable 
at all points, where x 2 (s) / xi(s). We claim that x 2 (s) / xi(s) for any s £ S (see Figure 4). Indeed, 
combining (66), (67) and (30), we get 

jx 2 (s) > tt > x 1 (s) for s £ [2k + 2, 2k + 3), 

\x 2 (s) < tt < xi(s) for s £ (2k + 1,2k + 2 - ^^], k £ N. 

Now consider the value s = 2k + 1, k £ N. In this case we have xi(s) = x 2 (s) = tt (see (31) and (55)). 
Notice that in this case we can not consider x 2 (s) as a root of the equation J(x, s) = (see Remark 
5.2 in Subsection 5.3). Therefore, let us return to the original function £b(x,s). Using the implicit 
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function theorem, we get x' 2 (s) = — 9 ^ 2 ^' s ^ f 9 ^*'^ \x=x 2 (s)- According to this formula, we calculate 
x' 2 {2k + l) for x = tt: 



x' 2 {2k + 1) = 

Note also that 
dg 2 {x,s) 



sec(A:7r) (-2 (kiT 2 (k + 1) + l) cos(2/ctt) + (2k + l)vr sin(2/cvr) +2) vr 



4A;(fc + l)(2fe + l)7rcos(A;7r) -4(3fc(A; + !) + !) sin(A:7r) 4£; + 2' 



d 



x 



sin(sx) — s((s 2 — l) xcos(x) + (s 2 — l) xcos(sx) + (s 2 + 3) sin(x) — 3ssin(sx)), 



^ 2 q^ S = 4 cos(x) + (— (s 2 — l) x 2 — 4) cos(sx) + x (—3 sin(x)s 2 + 2 sin(sx)s + sin(x)) 

are smooth functions and 

— -U =7 r = 8fe(l + fc)(l + 2A;)7r, ^ v ' ; | s=2fc +i = 4fc(l + fc)vr 2 . 

Thus X2(s) is continuously differentiable at s = 2k + 1. So we proved that 2:2(3) is a continuously 
differentiable function for Vs G S. Hence P2(m) is a continuously differentiable function for Vm G M. 

Continuity of the function 2:2(5) f° r s ^ (2A; + 2 — 2fc ^ 2 , 2/c + 2) follows from the implicit function 
theorem (see [13] p. 449). Indeed, the function J(x, s) is strictly monotone, thus 2:2(5) is a continuous 
function. 

Now we claim that there exist values s*, s e S such that X2(s*) = ir, 2:2(5) = xi(s), and s* < s. 
Indeed, 

x 2 (2k + 2 - —^2) < ^ < Xl ( 2A; + 2 - ^2) ' Xl( - 2k + 2) = ^ < a; 2 (2A: + 2), 

and xi(s) > ir Vs £ S. From continuity of X2(s) and xi(s), it follows that there exist s* and s as 
claimed. Therefore there exist the values m* and fh as claimed by Theorem 5.1. 

Finally note that at the point p = piijn) the function p \-¥ g2(p,m) changes its sign. Indeed, for 
m G (0, ^) this follows since the function G\(x) is increasing (see Subsection 5.1). For \ ^ m < 1 
this follows, for example, from the fact that g2(x,s) changes its sign over the interval (0, ^) (see 
Proposition 5.1). 

5.6. Study of P2{rn) for m > 1. We use the change of variables p = — , fh = — G (0, 1) to study 
Pi(m) in the case m > 1. We have g2(p,fh) = —^zg2(p, m )- Therefore, we have g2(p,fh) = iff 
92{p-,rn) = 0. In such a way we get the following functional equation: 

P2(m) = mp 2 (—). (69) 

Hence, the properties a) — e) (see Theorem 5.1) of p 2 {m) for m > 1 follow from the similar properties 
of P2(m) for m G (0, 1). 

6. Relative position of plots of p\(m) and P2(m) 

In this section we discuss the mutual behavior of the functions pi(m) and P2(m) and prove the 
following theorem: 

Theorem 6.1. The functions pi(m) and p2(m) have the following properties: 

a) pi(m) and p 2 ('m) are continuous functions for m G (0, 1) U (1, +00). 

b) For m G (0, 5) U (2, +00) there holds the inequality 

Pi(m) < p 2 (m). (70) 
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Fig. 4. Plots of the functions xi(s) and x%(s) (dashed line) 



For m G [i, 1) i/iere ftoZd i/ie following inequalities: 



pi(m) > pa(m) /or m G |±§ - 15+40fc+ % 2 ^ +8fc a 



Pi(m) < P2(Tn) for m G 



l+2fc l+fc \ 
3+2fc' 2+fc y 



Vfc G N. 



(71) 



P/o£s o/ t/ie functions pi(m) and P2(m) cross each other at the points m = -j^j and at points 
m = m k , where m k G (±±§ - 15+40fc+ 2 3 2fc^+8fc* ■ 3+11)' as Mows: 



Pi(m k ) = p 2 (m k ), pi 



k+l 



P2 



k + l 



irm 
1 — m 



■Kk, Vfc G N, 



For m G (1,2] there hold the following inequalities: 



Pi{m) < P2(m) for m G ( 
Pi(m) > P2{m) for m G 



( k+2 3+2fc 
fc+1' l+2fc 



3+2fc 



l+2fc + l+8fc(l+fc 2 ) ' fc y ' 



fc+1 



VA; G N. 



(72) 



Plots of the functions pi(m) and P2(m) cross each other at the points m = and at points 
m = m k , where fh k G (f$ff, f$H + 1+8 fc(i+fc^) ) ' as follows: 



Pi{m k ) = p2(m k ), pi 



k + l 



P2 



k + l 



k ) 1 \ k 
d) Both functions pi(m) and p2(m) tend to infinity as m 



7r(Jfe + l), V/cGN. 
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e) At the points m G {-^j, k G N} U {^T"l^ G N} l^ois of the functions pi(m) and P2{m) cross 
one another at an acute angle, since 

< ^(jfe + i 5 < Pl( ^TT } = + °°' ~°° = < P2{ — ] K °- 

Proof. Items a)-d) of Theorem 6.1 follow from Theorem 4.1 and Theorem 5.1. First notice that proof 
of the case m > 1 is reduced to proof of the case m G (0,1), since the functions pi(m) and P2(jn) 
satisfy the functional equations pi(m) = rnpi(^) and P2{rn) = mp2(^). Further, if m G (0, |), then 
Pi(m) < 3™ < 5.7m < p 2 (m). If m € [3, §), then Pl (m) ^ ^ < P2M (see (51)). If m G [±, 1) 
then we have the following inequalities: 

P2M < ^ < pi(ro) for m G - ig+aoF^F+IF 

P2(m)>S>PiHforme [^f,^), ^N. 
Since pi(m) and P2(m) are continuous functions, then these inequalities imply existence of values 

_ _ /1 + 2/c 2 1 + 2F 

mfc G V3 + 2/c ~ 15 + 40/c + 32k 2 + 8/c 3 ' 3 + 2/c / 
such that p\(fhk) = P2{rnk) for any G N. It follows immediately from Theorems 4.1 and 5.1 that 

Pi(iTT)=P2(iTT) = S, VfcGN. 

Finally, item e) is checked by direct calculations. This completes the proof of Theorem 6.1. □ 

Figure 5 shows the plots of the functions pi(m) and p2(m). Note that these plots have an infinite 
number of intersection points. 

7. Limit behavior of Maxwell sets and cut time 

Consider a sequence of extremal trajectories Qt as p — > 0, m — > in. The instants of time t = 
t\{m) = 2pi(m)/m and t = t2(m) = 2p2(m)/m define asymptotics of the first Maxwell time for Q t , 
corresponding to the reflections e 1 and e 2 of mathematical pendulum (6). We obtained two-sided 
estimates of ti(m) and t2(m) for m G (0, 1) U (1, +00) and discussed continuity and differentiability of 
these functions. Note that the functions t\{m) and t2(m) are characterized by very complex behavior: 
their plots have vertical tangents at a countable number of points and lim t\(m) = lim t2(m) = +00. 

Then we showed that plots of t\{rn) and t2(m) have an infinite number of intersection points. 

Denote by A the vector of adjoint variables. In the paper [7] limit behavior of the Maxwell set 
MAX 1 is studied. Then upper bound on cut time 

^cut(A) = sup{t > I Q s = Exp(A, s) is optimal for s G [0, t]} 

as (9, d) —7- (0, 0) was obtained. Namely there is the following estimate for the cut time: 



lim tcut(A) ^ ti(m) for m > 0, m / 1. 

p— ¥0,m— >m 

Fix any compact set K C {m G M. \ m > 0, m/1}. In [7] the following upper bound for cut time is 
proved: 

lim t cut (X) ^ maxii(m). 

p—t0,m£K m&K 

Similar estimates for the cut time i cu t, based on the limiting behavior of the Maxwell set MAX 2 , will 
be obtained in later studies. 

Recall that plots of t\{m) and t2{m) have an infinite number of intersection points. This means 
that even in asymptotic case the behavior of the first Maxwell times for the plate-ball problem is much 
more complicated than for the related invariant optimal control problems (nilpotent sub-Riemannian 
problem with the growth vector (2,3,5) [15], the problem on Euler elasticae [16], sub-Riemannian 
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12 3 4 



Fig. 5. Plots of the functions p\{m) and P2(m) (dashed line) 

problem on the group of motions of a plane [17], sub-Riemannian problem in Martinet case [14]). The 
first three problems were studied by Yu. Sachkov. He showed that in these problems the plots of 
Maxwell times have not more than two intersection points globally. While in the plate-ball problem 
this number is infinite even in the simple asymptotic case. This shows the new level of complexity 
of the plate-ball problem not encountered in control theory before. Taking into account complexity 
of parameterization of extremal trajectories in this problem, it is difficult to obtain the exact solu- 
tion. However, on the basis of these results it is possible to develop algorithms and software for the 
approximate solution of the plate-ball problem. This will be the subject of future work. 
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